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ABSTRACT 

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey 
(SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the 
expected number of quasars lensed by early-type galaxies and their image separation distribution 
assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 
to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant 
model {w = —1) we obtain J^a — 0.741q :[g(stat.)lg'|!)g(syst.). Allowing w to be a free parameter 



we find ft 



M 



0.26t5];!]^(stat. 



'^ao^syst-) and w 



-1.1 ± 0.6(stat.)lQ'5(syst.) when combined with 



the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy 
sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, 
and provide additional confirmation of the presence of dark energy consistent with a cosmological 
constant, derived independently of type la supernovae. 

Subject headings: cosmological parameters — cosmology: theory — gravitational lensing 
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1. INTRODUCTION 

The accelerating expansion of the universe is one of 
the central problems in modern cosmology. This accel- 
eration is usually attributed to the dominant presence 
of a negative-pressure component that is often referred 
to as dark energy. There are many models that explain 
the acceleration, including a classical cosmological con- 
sta nt, decaying scalar field s, and topological defects (e.g., 
see [Peebles fc Ratral 120031 for a review). In addition it 
might also be explained by lon g-range modificatio ns of 
the gravitational force law (e.g.. ICarroU et al]l2004f ). 

The dark energy is characterized by its cosmological 
density f^DE, and its equation of state w, which is de- 
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fined as the pressure divided by the density of dark en- 
ergy. In particular, measuring w is a useful test of mod- 
els for dark energy. Since w determines the expansion 
rate of the universe and the cosmological distance to a 
given redshift, not only Que but also the value of w 
and its time dependence can be inferred from distance 
(or volume) measurements on cosmological scales. One 
powerful way to constrain w is via observations of dis- 
tant type-la supernovae that make use of luminosity- 
decli ne rate correlations to standardize their luminosi- 
ties (|Riess et al.l Il998t iPerlmutter et al.l [1999' ): Since 
the standardized luminosities of type-la supernovae have 
small scatter, they serve as an excellent standard can- 
dle to measure cosmological distances. Another probe 
of dark energy is the fluctuation sp ectrum of the cos- 
mic r nicrowave background (C MB; Ide Bernardis et aTl 
I2000t ISpergel et al.|[2003L 120071 ). The integrated Sachs- 
Wolfe effect, which can be detected by correlating the 
CMB map with the large-scale distribution of galax- 
ies, allows a direct detection o f the dark energy com- 
ponent (e.g.. iRassat et al]|2007D . Additional constraints 
on dark energy come from bary on acoustic oscillations 
in the galaxy power spectrum (lEisenstein et al.l l2005l : 
ICole et al.| [ 2005l: iPerciyal et al.ll200'jD and X-rav c lusters 



of galaxies (|Allen et al.ll2004 120071: iRapetti et al.ir2005D . 

Since different methods involve different systematics and 
degeneracies with the cosmological parameters, it is of 
great importance to use as many independent observa- 
tions as possible in studying dark energy. 

The statistics of strong lensing offer an alterna te con- 
strain t on dark e nergy ([T urner 1990; Fukugit a et al.l 
119901 but see also lKeetoni r2002) . The probabihty that 
a distant object is strongly lensed is proportional to the 
number of possible lensing objects along the line of sight, 
and thus quite sensitive to dark energy. This method 
has been applied to both optical and radio lens samples 
to derive interesting constraints on the value of the cos- 
mological constant (iMaoz fc Rixl 1993 ; iKochanefl 119961 : 
iFalco et al.lll998l : IChiba fc Yoshii|[T999l ). but such appli- 
cations have been limited by the small size of existing 
lens samples as well as poor knowledge of source and 
lens populations (|Maod [20051 ). For instance, past work 
tended to rule out large values (> 0.7) of Qbe (e-g-, 
lKochane^ll996l ). because of ove restimates of the lum i- 
nosity function of galaxies (e.g., IChiba fc Yoshii|[l999D . 
The most recent lensed quasar survey in th e radio band, 
the Cosmic Lens All-Sky Survey (CLASS; iMvers et all 
l2003t iBrowne et al.ll2003[ ). contains a statistical sample 
of 13 lenses. Cosmological constraints from this lens 
sample are roughly consistent with the current standard 
model in which the universe is dorninated by dark energ y 



(|Chae et al.ll2"00l IChai[200a [2007t iMitchell et alll200l 
The Sloan Digital Sky Survey Quasar Lens Search 
(SQLS;[0 guri et all2006( l provides a large statistical lens 
sample appropriate for studying cosmology. It is based 
on the optical quasar catalog from the Sloan Digital Sky 
burvey fSDSS: lYork et al.l[2000f l. and therefore is comple- 
mentary to the CLASS in several ways. In particular the 
well-known redshift dis tribution of quasar s and the lens- 
ing selection function (jOguri et al.l |2006| . hereafter Pa- 
per I) allow an accurate estimate of lensing rates. We 
present our fir st complet e lens sample from Data Release 
3 fDR3: lAbazaiian et al. 2005: Schneider et aL 2005) in 
llnada et al.l ( 2007L hereafter Paper II): It consists of 11 



lensed quasars with flux ratios of faint to bright images 
greater than 10^°-^ (for double lenses) and image sep- 
arations between 1" and 20", selected from 22,683 low- 
redshift (0.6 < z < 2.2) quasars brighter than i — 19.1. 
In this paper, we use a subsample of this optical lens 
sample to constrain cosmological parameters, in partic- 
ular the dark energy abundance and equation of state. 

This paper is organized as follows. In ij2| we briefly 
summarize our statistical lens sample. Section [3| de- 
scribes how the expected number of lensed quasars is 
computed. We present our results in ijH and summa- 
rize in Sj5l We denote the present matter density as 
Qm- The present dark energy density is described as 
r^DE, or Oa if the cosmological constant w = — 1 is as- 
sumed. We use the Hubble constant in dimensionless 
form h — i/o/(100kms~^Mpc~^). Throughout the paper 
we assume a flat universe Qm + Qde — 1- Magnitudes 
quoted in the paper are corrected for Galactic extinction 
(jSchlegel et al.l[l998h . 

2. LENSED QUASAR SAMPLE 

Our lensed quasar sample is constr ucted from the 
SDSS DR3 spectroscopic quasar catalog ([Schneider et al.l 
l2005f l. The properties of the SDSS are pres ented in a 
series of technical papers. iGunn et al.l (|2006f) describes 
the dedicated wide-field 2.5-m t elescope. Deta ils of the 
photometric survey a,re given in iFukugita et a l. (1996), 
iGunn et all (Il998l). iLupton et al.l (Il999l). iHogg et all 



(20011). iLupton et all (|2001D. ISmith et al.l (I2002D. 
lPieretal.1 f 20031 'Ivezic et al.' ('2004'), 'Tucker^etalJ 
(|2006l ). and Lupton (2007). Blanton et al. (2003) present 
the tiling algorithm of the spectroscopic survey. Spectro- 
scopic quasar tar gets are selected accor ding to an algo- 
rithm described in lRichards et al.) (|2002| ). Details of each 
public data set arc given in a scries of data release pa- 
pers (Stoughton et al. 200 2; Ab azajian et al. 2003, 200|, 
[20051 : rAdelman-McCarthvll2006L [2007i ) . 

The DR3 statistical lens sample contains 11 lensed 
quasars (see Paper II) . The sample is restricted to a range 
of i-band flux ratios, which are the fluxes of the fainter 
images divided by those of the brighter images for dou- 
ble lenses (no condition on flux ratios is set for quadruple 
lenses), fi > 10~°-^ and image separations 1" < 9 < 20" 
where the completeness of the candidate selection is al- 
most unity (see Paper I). In this paper, we apply two 
additional cuts to select a subsample appropriate for our 
dark energy study. First, we restrict the image separa- 
tion range to 1" < 9 < 3". At 9 < 3" lens potentials are 
in many cases dominated by those from individual lens- 
ing galaxies, whereas beyond 9 = 3" the contribution 
of surrounding dark matter t o lens potentials b egins to 
become more significant (e.g., iKochanek fc Whitd 120011 : 
lOguri et al.l[2005bl : [Ogurl [2006|) ■ The effect of the ex- 
ternal field can, in principle, be included in our theo- 
retical model of lensing rates. However, it is difficult to 
observationally constrain the probability distribution of 
external fields (note that we adopt observationally de- 
termined velocity functions for our computations. See 
§3.ip . indicating that it introduces additional systemat- 
ics. Second, we require that the lensing galaxy be fainter 
than the quasar components. If the lens galaxy is too 
bright, it will strongly affect the colors of the quasars and 
the lensed quasar wi ll not be targeted for spectroscopy 
([Richards et al.|[200^ . biasing against us discovering the 
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TABLE 1 
SDSS DR3 Quasar Lens sample 



Name 


A^img 




2; 


^cor 


e 


Lens 


Note 




Ref. 


SDSS J0246-0825 


2 


1.685 


0.723 


17.77 


1.04 


E 






1, 2 


SBS0909+523 


2 


1.377 


0.83 


16.17 


1.11 


E 






3, 4 


SDSS J0924+0219 


4 


1.523 


0.393 


18.12 


1.78 


E 






5, 6, 7 


SDSS JlOOl+5027 


2 


1.839 




17.31 


2.86 


E? 






8 


SDSS J1021+4913 


2 


1.720 




18.97 


1.14 


7 






9 


SDSS J1226-0006 


2 


1.125 


0.517 


18.23 


1.24 


E 






10, 11 


SDSS J1335+0118 


2 


1.571 


0.440 


17.53 


1.57 


E 






12, 11 


Q0957+561 


2 


1.413 


0.36 


16.67 


6.17 


E 


e > 3" 




13, 14 


SDSS J1004+4112 


5 


1.740 


0.68 


18.84 


14.6 


C 


> 3" 




15, 16 


SDSS J1332+0347 


2 


1.438 


0.191 


17.89 


1.14 


E 




> 


17 


SDSS J1524+4409 


2 


1.210 


0.310 


18.76 


1.67 


E 


^qso ~ "igal 


> 


18 



References. — ( 1) llnada et aLll2005l: ( 2) lEiKcnbrod ct al. 200 71: f3l lOscoz et aL|[T997| : 
(4)'Lubin et al."2000'; (5)'Inada et al."2003a'; (6) Ofek ct aL 2006; (7VEigenbrod et al.'200'6al: 
^8) Osuri et al. 2005a; (9) Pindor et al. 2006; (10) Inada ct al. 2003c; (11) EiKcnbrod c t all 
2006ta; (12) iQguri et al n2004b; (13) iWalsh et al.l |1979|: (14) lYoung et al.nl981j; (15) 



Inada et al ] |2003bl : fl6Y lOguri et"ani2004al : (171 IMorokuma et al.ll2007l: fl8) lOguri et al.l 
20071 



Note. — See Paper II for details of the construction of the statistical lens sample, icor 
is the i-band PSF magnitude of the object with SDSS spectroscopy, corrected for Galactic 
extinction. The image separation 8 is defined by the maximum separation between any 
image pairs. "Lens" indicates the morphology of the lensing galaxy: S=spiral; E=elliptical; 
C=cluster; ?=unknown. We adopt the first 7 lenses for our statistical study; the 4 lower 
lenses are not used because of reasons indicated in the Note. 



lens. In addition, lens galaxy fluxes add to the bright- 
ness of the system, which could enhance the number of 
lenses in the flux-limited sample. These biases make the- 
oretical predictions much more difhcult and uncertain. 
Therefore, we require that the i-band point spread func- 
tion (PSF) magnitude for the quasar components, Zqso, 
be brighter than the i-band magnitude of the lensing 
galaxy, Zgai. These two cuts remove four lensed quasars 
from the lens sample. Hereafter, we use the remaining 
seven lensed quasars to constrain cosmological parame- 
ters. Table [T] summarizes the lens sample we use in this 
paper. 

3. THEORETICAL MODEL 

In this section, we describe how to compute the ex- 
pected number of lensed quasars in the SQLS DR3 
sampl e in a cosmological model, following [Turner et al.l 
with the selection function taken into account. 

3.1. Lens Galaxy Population 

We consider early-type galaxies as lensing objects. 
Although late- type galaxies are more abundant, stan- 
dard models predict that the strong lensing probabil- 
ity is domin ated by that from e arly-type galaxies (e.g. , 
Turner et a l. 1984; iMaoz fc Ri^ Il993t iKochanekl 119961 : 



MoUer et al...2007. ). This is particularly true if we restrict 
image separations to be larger than 1", because late-type 
galaxies have smaller velocity dispersions on average and 
therefore have smaller mean image separations. Indeed 
this is confirmed by observations. Only a few of the > 60 
known lensed quasars with 9 > 1" are produced by late- 
type galaxies. "^^ Moreover, none of the lensed quasars 
in our sample appear to be caused by late-type galaxies 
(see Table II]). 

The contribution of large-scale dark matter fluctu- 
ations around lensing galaxies (environmental conver- 



33 gee the CAST LES 

|http: / / cfa-www.harvard.edu/castles/ 1 



webpage 



at 



gence) is important because i t could bias cosmological re - 
sults from lensing statistics (|Keeton fc Zabludofil 12004'). 
The maximum image separation of 3", however, makes 
the effect of external convergence due to associated dark 
matter on the lensing probability moderate. Specifically, 
the environmental convergence enhances th e integrated 
lensing probability only by < 10% at 6* < 3" (jOguri et al.l 
i2005b,) . Therefore we neglect external convergence, al- 
though we examine its impact on our results below. 

It has been argued that the radial mass profile of 
galaxies can be approximated as a singular isothermal 
sphere for the scales re levant for strong lensi ng (e.g., 
iRusin fc Kochanek 2005tlKoopmans et al.ll200"6l ). In this 
paper, we adopt an elliptical version of this, a Sin- 
gular Isothermal Ellipsoid (SIE). The ellipticity does 
not significantly aff ect the total lensing cross section 
(jHuterer et al.|[2005h , but including ellipticities allow one 
to take account of the different selection functions of 
double and quadruple lenses (see Paper I). The two- 
dimensional surface mass density of an SIE with elliptic- 
ity e at a position x and y from the center of the galaxy, 
with the X axis aligned with the major axis, is given by 



k(x,?/) 



=A(e) 



1-e 



(1 — e)^x-^ + y 



1/2 



(1) 



where denotes the Einstein radius, which in turn is 
related to the galaxy velocity dispersion cr^ by 



0^ = in 



V C / Dns 



(2) 



where Dig and Dos are the angular diameter distances 
from lens to source and from observer to source, respec- 
tively. The parameter A(e) is the velocity dispersion nor- 
malization factor for non-spherical galaxies. 

It is not straightforward to determine the normaliza- 
tion factor A(e). What is needed is the calculation of 
velocity dispersions for lens galaxies, which depends on 
three-dimensional shape of the lensing galaxies when the 
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assumption of spherical symmetry is relaxed. Two ex- 
treme possibilities are that all galaxies have either oblate 
or prolate shapes. In this paper, we assume that there 
are equal number of oblate and prolate galaxies and com- 
pute the normalization by taking the avera ge of norma l- 
izations in the oblate and prolate cases (seelChaillOOl). 
This is a reasonable assumption given the fact that un- 
derly ing dark halos have triaxial shapes (|Jing &: Sutol 
I2OO2I) . Moreover the observed axis ratio distribution is 
con sistent with a po pulation of triaxial early-type galax- 
ies (jVincent &: Rvden 2005). For the distribution of el- 
lipticities, we adopt a Gaussian distribution with mean 
e = 0.3 and dispersion de — 0.16 (but truncated at 
e = and 0.9) that is consistent with the observed el- 
lipticity distri bution s of the light of early-type galaxies 
(iBender et al.l 1X9891: iSaglia et al.lll993t [Jorgensen et al 



1995 : iRest et al. l l200lHAlam fc Rvdenll2002l : ISheth eta 



20031) . 

One of the most important elements in predicting the 
number of lensed quasars is the velocity function of galax- 
ies. As a fiducial velocity functi on we adopt that derived 
from the latest SDSS DR3 data (|Choi et al.ll2007D . which 
is fitted by a modified Schechter function of the form 



dn 
day 



exp 



P dUy 

r(a//3)"^7' 



(3) 



where ct*, a, /3)=(8.0 x IQ-^h^M^c'^ , 161kms-\ 
2.32, 2.67). This is some what different from the ve- 
locity function deri ved by ISheth et"al1 ()2003l see also 
IMitchell et all l2005l) upon SDSS DRl data. They fit 
the same functional form, but found different parame- 
ters, ((/)*, CT*, a, /3)=(4.1 X 102iMM2C"^88.8kms-\ 
6.5, 1.93). We will use the ISheth etall (|2003D parame- 
ters to estimate the size of the systematic error due to 
the velocity function. 

Since the lens galaxy will typically have D\s ~ 0.5Dos, 
the lenses in our survey can have redshifts up to z ~ 
1 (e.g., SBS0909-K523 has the lens redshift of 0.83), 
and thus any redshift evolution of the velocity func- 
tion, which could chan ge the lensing rate and degen- 
erate with cosmology ()Keetonl |2002|) . must be taken 
into account. While it has been argued that early- type 
galaxies evolve onl y through passive luminosity evolu- 
tion at z < 1 fe.g..' llm et aLll2002t Iwillmer et al.ll2006f ). 
theoretical studies fav or slight evolution with redshift 
through mergers fe.g.. 'New man fc Davis! 12000 ^). On the 
other hand the lens redshift distribution of strong lens- 
ing i s consistent with no evolution of the velocity func- 
tion (I Of ek et al.l[200llChae fc Maol [200l IMitchell "eTaH 



l2005t iCapelo fc Nataraianll2007t l. In this paper we as- 
sume that the velocity function does not e volve, b ut we 
also consider the evolution model used in [Chad (120071 
based on a semi-analytic model bv lKang et al.ll2005[ ) as 
well to estimate the systematic impact of the evolution 
on our result. In the model the number density and the 
characteristic velocity dispersion are simply replaced by 
0* 0*(1 + 2;)-0-229 and cr, ^ cr,(l -|- z) 



-0.01 



3.2. Quasar Luminosity Function 

The quasar luminosity function is used to calculate 
the magnification bias. We adopt the luminosity func- 
tion constrained from the combinat ion of the SP S S and 
2dF (2SLAQ), namely the 2SLAQ-f jCroom et all (pOOl ) 



model in lRichards et al.l ()2005[ ). as our fiducial model: 



$(Afg) 



lQaA(l-M(M,-M;) _^ ^Q0.4(l-A)(Mg-M;) ' 



M;{z) = M*(Q) - 2.5(fciz + k2Z^), 



(4) 
(5) 



with the parameters of (/3h, A) -^g(0)i 
fc2) = (3.31, 1.45, 1.83 X 10-^(/i/0.7)3Mpc-Vag-i, 
-21.61 + 51og(/i/0.7), 1.39, -0.29). The bright and 
faint end slopes are broadly consis tent with those in the 
bolometric luminosity function of iHopkins et al.l (120071 ) 
at z ~ 1 — 2. The luminosity function is in terms of 
rest-frame 5-band absolute magnitudes: We convert it 
to observed i-band a pparent magnitudes u sing the K- 
correction derived in iRichards et all ()2006t ). Since the 
luminosity function was derived assuming = 0.3 and 
^Ia — 0.7, we adopt this cosmology for computing the 
absolute magnitudes used to compute the magnification 
bias no matter what values of fJjvi and w we consider for 
the remainder of the analysis. 

3.3. Number of Lensed Quasars 

The lensing cross section CTicns for a given lens is 
computed numerically using the public code lensmodel 
P<:eetonll200ll) . We compute the sum: 



= / du 



/ia>(i) ' 



(6) 



over the source plane positions u (with magnification /x) 
where multiple images are produced. The cross sections 
are computed in units of 0e and they are weighted by the 
ratio of the differential luminosity functions in order to 
take magnification bias into account. The suffix i indi- 
cates the number of images, with i = 2 for double lenses 
and j = 4 for quadruple lenses. For double lenses the 
integral is performed over the region where the flux ratio 
of faint to bright images is larger than 10~°'^, in order to 
match the selection function of our lensed quasar sam- 
ple. In Paper I, we found that the magnification factor 
of lensed images depends on the image separation, be- 
cause the SDSS quasar target selection adopts the PSF 
magnitude for the magnitude limit. It was found that 
the dependence is fitted well by the following form (see 
Paper I for details): 



and 



= mot + (1 - M)Mbri, 



^ = - [1 + tanh(1.76 - 1.786*)] 



(7) 



(8) 



where 9 is in units of arcsecond and /itot and /ibri are 
the total magnification and magnification of the bright- 
est image, respectively.'^^ We compute the lensing cross 
section as a function of dimensionless image separation 
= 6'/6'e, i.e., daicns,i/d9. 

From the lensing cross section, we can compute the 
differential probability that a source ed z — Zs is strongly 

•^■^ We note that eq. (15) of Paper I holds only approximately if 
the slope of the source luminosity function is close to —2. In this 
paper we compute the full magnification bias using eqs. (|6j and 
0. 
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lensed with the image separation 6 as 



dpi 

de 



s ; f'qso J 



I 



cdt 
dzi——[l 
dzi 



dn 

day- — 



de{D,,9^) 



de 



dO dffy dn 
§ d0E day 



, cdt^ 
dzi——{l 

dzi 



[DM' 



I daic 



dO 



sal 



so), (9) 



where 4ii{9) is the conipleteness of our lens candidate se- 
lection estimated from simulations of the SDSS images 
(see Paper I; for double lenses we adopt the complete- 
ness averaged over the flux ratio between 10~°'^ and 1), 
cdt/dzi denotes the proper differential distance at zj, 
and (-DoI^e)^ converts the lensing cross section from the 
dimensionless unit to the physical unit. The Heaviside 
step function 0(igai — iqso) is added to include the con- 
dition that the quasar components should be brighter 
than the lensing galaxy. We compute igai using the 
correlation between the luminosity and velocity disper- 



sion o f early-type galaxies measured by Bcrnard i et al 



2003) combined with K-correction from Fu kugita et al 



..1995f). The effect o f the luminosity evolution measured 
bv iBernardi et ahl (|2003f ) is included. The estimated 
galaxy magnitudes are broadly consistent with those ex- 
pected from the empiric al scaling observed for known 
lenses (jRusin et al.ll2003[ ). The total lensing probability 
for image separations between 1" and 3" is then given by 



s ; i-qso J 



dO—[Zs,lqso)- 



(10) 



The expected number of lensed quasars in our quasar 
sample is computed by counting the number of quasars 
weighted by the lensing probability. Ultimately this can 
be done by adding the probabilities for all source quasars: 

Ni ^ ^ Pi{Zs,iqso)- (11) 
source QSOs 

To save computational time, we actually calculate the 
expected number of lensed quasars for each redshift- 
magnitude bin and then sum over the bins. If 
A'qso(^sj , *qso,fc) IS the uumbcr of source quasars in the 



redshift range 
magnitude range 



Azs/2 < Zs < Zsj + Azs/2 and a 
Ai/2 < iqso < iqso,fc + Ai/2, 



then the predicted number of lensed quasars is 



A^qso (^s j' 7 '^qso, k^Pi (^s,j : ^qso.fc) • (-^^) 



We adopt bin widths of Azg — 0.1 and Ai — 0.2. 

3.4. Likelihood 

We perform a likelihood analysis on the DR3 SQLS lens 
sample. The likelihood i s corn puted using the method 
introduced bv iKochanekl (|1993l ) 



In L = In 



lens 



lens unlensed QSOs 

dpi 



de 



{N2 + Ni), 



(13) 




Fig. 1. — Relative likelihoods of the value of the cosmological 
constant Ha from fitting the SQLS DR3 data, assuming a spa- 
tially fiat universe. The vertical dotted lines indicate the Icr range 
estimated from Ax^ = 1. The likelihood becomes maximum at 
= 0.74. 

where dpi/de is calculated from eq. ^ and N2 (dou- 
bles) and N4 (quadruples) are from eq. (fT^ . We note 
that the valid approximation P27 P4 ^ 1 is used here. 
We neglect three image events caused by naked cusps 
because they make a negligible contribution to the total 
cross section. The summation in the first term runs over 
the seven lensed quasars listed in Table [T] The distribu- 
tion of lens redshifts offers an independent test of cosmo- 
logical model, however it has been shown that it is more 
sensitive to the redshift evolution of the lens galaxy pop- 
ulation as well as the select ion bias (e.g. JOfek et al.ll2003l : 
ICapelo fc Nataraianll2007l ). We do not include informa- 
tion on the lens redshift zi (i.e., we adopt the probability 
after integrating over the lens redshift in eq. [5] ) because 
the lens redshifts are not known for all lenses (see Table 
[1]) and including zi for those objects that do have lens 
redshifts might introduce systematic effect related to the 
incompleteness of our redshift information. 

4. RESULTS 

In this section, we derive constraints on cosmological 
parameters from the lens statistics. Since the abundance 
of lensing galaxies we adopt in this paper was determined 
observationally, our results do not depend on the normal- 
ization (Tg or shape of the primordial power spectrum. 
The Hubble constant is not important for these calcula- 
tion because no absolute length scale is used. We assume 
a flat universe flM + ^de = 1, therefore there are only 
two independent cosmological parameters in our analy- 
sis, fluE {^a) and w. 

4.1. Cosmological constant 

First we derive constraints on the cosmological con- 
stant r^A assuming w = —1. We compute the likelihood 
(eq. [13]) as a function of f^A, which is plotted in Figure 
[TJ The resulting constraint SIa = 0.741° ^5 (Icr; the error 
is estimated from Ax^ = — 21n(L/Lmax) — 1) is broadly 
consistent with oth er measurements o f the cosmolog- 
ical constant (e.g.. ISpergel et all l2007l and references 
therein). It is also consistent with CLASS, for whi ch the 
best-fit cosmological constant w as l^A = 0.7-0.8 (ICha^ 
[200a \200% IMitchell et all 120051) . The case J^a = has 
Ax^ ~ 9 (the expected total number of lenses in the DR3 
sample is ^ 1) and is therefore rejected at the 3a level. 




^T-t-HH 

1 2 3 

S[arcsec] 

Fig. 2. — The number distribution of lensed quasars is plotted 
as a function of the image separation 6. The histogram shows the 
number distribution in the SQLS DR3 statistical sample (see Table 
[T] The bin size is O'.'S, thus the actual number of lenses in each bin 
is half of what we plot). We use only lenses in the image separation 
range 1" < < 3" as indicated by the vertical dotted lines. The 
solid line indicates the prediction of our best-fit model Q\ = 0.74 
(see Figure [T]l. The dashed line shows the predi ction of our b est-fit 
model when we adopt the velocity functio n of ISheth et al.l (120031) 
i nste ad of our fiducial velocity function of lChoi et al.! (2007). See 
ii4.3l for a detailed discussion of the effect of adopting the different 
velocity functions. The sharp decline below 9 = 1" is due to the 
selection function 4>i{9), which rapidly decreases at 6 < 1". Note 
that our statistical lens sample contains two more lensed quasars 
at > 3" that are not shown in this figure. 

We test the validity of our best-fit model (J^a — 0.74) 
by comparing the image separation distribution with the 
observation. In Figure [5] we plot the expected number 
distribution (eq. [9] summed over all sources) with the 
image separations of 7 lensed quasars in the statistical 
sample. Although the small number of lensed quasars 
prevents detailed comparison, both the normalization 
and the overall shape of the curve appear to match the 
observed one. The Kolmogorov-Smirnov test finds that 
the best-fit model is consistent with the observed distri- 
bution at a significance level of 72%. 

4.2. Dark energy 

Next we derive constraints on the equation of state w 
as well as the dark energy abundance r^DE- To do so we 
compute the likelihood as a function of Om = 1 — f^DE 
and ly, and compute constraints in the two-parameter 
space. The result is shown in Figure [3l The confi- 
dence region shows degeneracy in a similar direction as 
constraints from type la supe rnovae (e.g., lAstier et al.l 
l2006t IWood-Vasev et aLll2007[ ). therefore strong lensing 
alone does not constrain these parameters very well. As 
a complementary constraint, we also consider the likeli- 
hood from the measurement of the scale of baryon acous- 
tic oscillations ( BAO) in the SPSS lum inous red galaxy 
power spectrum (jEisenstein et al.ll2005( l. Specifically we 
adopt their constraint A = Dy (0.35) y/ ft mHq /0.35c = 
0.469 ± 0.017, where £)y(0.35) is the dilation scale at 
z — 0.35. From the combined constraint shown in Fig- 
ure H we obtain ilu = 0.26tom and w = -1.1 ± 0.6 
(Icr). Again the constraint is c onsistent with other 
measureme nts of dark energy (e.g.. Tegmark et al .ll2006l: 
Astier et al...2006. : .Spcrgcl et al...2007, : .Wood-Vasev et al.l 




0.2 0.4 0.6 



Fig. 3. — Contours at la and 2cr confidence levels (estimated from 
Ax^ = 2.3 and 6.17) are plotted in the Hm-w plane. Solid lines 
indicate the constraint from the SQLS DR3, whereas dotted lines 
are from the baryon acoustic oscillations (BAO) detected in the 
SDSS luminous red galaxy power spectrum (Eisenstein et al. 2005). 
The joint constraint from SQLS and BAO arc shown by shaded 
regions: The best-fit model (Qm, ui)=(0.26, —1.1) is indicated with 
a cross. 

4.3. Systematic errors 

In this paper we have made a number of assumptions to 
compute the expected number of lensed quasars in each 
cosmological model. Here we examine the sensitivity of 
our results to these assumptions. To do so, we change 
the input model within the expected uncertainties, and 
compute the best-fit cosmological parameters to deter- 
mine their sensitivity to these details. We consider the 
following sources of systematic errors. 

• For the mass distribution of the lens galaxy, we 
change the fraction of prolate/oblate shapes by 
±25% to derive a rough estimate of the systematic 
error coming from the dynamical normalization. 
We also change the peak of the ellipticity distribu- 
tion by ±0.1, roughly correspo nding to the current 
uncertainty in o bservations (Bender et al. 19891 
Saglia et al.lll993l:Dorgenseii et al. 1995- Rest et al] 
20011 : lAlam fc Rvdenll2002HSheth et al.ll2003D . 

• We change the faint end slope of the quasar lumi- 
nosity function by ±0.2 while keeping the other pa- 
rameters fixed, in order to examine the sensitivity 
to the quasar luminosity function. We experiment 
with the faint end slope because the most poorly 
constrained part of the quasar luminosity function 
is the faint end slope with sort ie evidence that it 
is shallower than in our model () Jiang et al.| [20061. 
The observed quadruple lens fraction in the SQLS 
also f avors slightly shallower faint end slope (,Oguril 
12001 . 

• In this paper we have neglected the contribution 
of external convergence and shear (jOguri et al.l 
|2005b): To see how the external fields change the 
predicted number of lenses we simply assume the 
redshift-indepcndent distributions of external fields 
are 0.05 ± 0.2 dex for shear and 0.01 ± 0.5 dex for 
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TABLE 2 

Systematic errors 



Uncertainty 


w = —1 


wjt-l (with BAO) 


ASJa 




iXw 


prolate 25% - 75% 


±0.04 


±0.02 


+0.1 
-0.2 


e ±0.1 


-0.01 


+0.00 


+0.0 


01 ^ ±0.2 


±0.04 


±0.02 


±0.2 


external shear 


+0.00 


-0.00 


-0.0 


external convergence 


-0.01 


+0.01 


+0.1 


Sheth et al. dn/da^ 


+0.10 


-0.04 


-0.4 


dn/doy evolution 


+0.04 


-0.02 


-0.2 


^qso ^gal ^ lbO.5 


+0.04 
-0.02 


±0.01 


±0.1 


total 


+U.i3 

-0.06 


+U.U3 
-0.05 


+0.3 
-0.5 



Note. — Total errors are estimated from the quadrature 
sum of all errors. 



converge nce (e.g iDalal fc Watson|[200l see also 
rMomchc vaet al.l |2006[) . 



• For the velocity fun ction, we rep l ace th e function 
we use with that of ISheth etall (|2003l ). We also 
consider the effect of redshift evolution of the ve- 
locity function by using a simple parametric model 
(see §3.11 for details). 



The condition 



^qso 



< may involve some 



uncertainties in estimating the galaxy luminosity or 
the condition itself. Thus we allow an uncertainty 
of ±0.5 in the cut to estimate the systematic error 
that this introduces. The lens sample is unaffected 
by this change. 

Our results summarized in Table [2] indicate that the 
largest uncertainties in our conclusions come from the 
dynamical normalization, the faint end of the quasar 
luminosity functions, the velocity function of the lens 
galaxies and its redshift evolution. In particular, chang- 
ing the velocity function significantly increases the best- 
fit va l ue of the cosmological constant, as discussed by 
[CI^(|2n0l . In essence, the lensing optical depth scales 
as (7^, so small errors in the velocity function lead to 
much larger errors in the cosmological estimates. Ide- 
ally we would self-calibrate t he velocity funct ion based 
on the image separations (see lKochanek|[r993() . but the 
two models we consider here have very similar predic- 
tions for the image separations (see Figure [21) and the 
data only favor the Choi ct al. (2007) model by the 
difference of ~ 0.5. Our result indicates that the errors 
from the galaxy ellipticity and external perturbations are 
negligibly small. If we combine all these uncertainties, 
we obtain systematic errors comparable to the statisti- 
cal errors, e.g., 0.^ = 0.74lo it5(stat.)lo |!)g(syst.) for the 
cosmological constant case. 

5. SUMMARY AND DISCUSSION 

We have derived constraints on dark energy using the 
new optical strong lens sample from the SQLS DR3. 
We take various selection effects into account to make 
reasonably robust predictions for the number of lensed 
quasars in the SQLS. We have found that the de- 
rived constraints agree well with the current concordance 
cosmology. Assuming a cosmological constant (w — 



-1) we have obtained ^Ia = 0.74_Q'j^5(stat 



)^o;o6(ByBt. 



The constraint primarily arises from the total num- 
ber of lenses in our statistical sample. For the more 
general case w ^ —1, the constraint was combi ned 
with that of the SDSS BAO (lEisenstein et al.ll2005f) to 
break the degeneracy. The resulting joint constraint 
is = 0.26lJ^;^^(stat.)l°;°^(syst.) and w = -1.1 ± 
0.6(stat.)t^J(syst.). The results are in good agreement 
with recent con straints from radio lenses (Chac 2003|, 
120071: [Mitchell et al . 2005). The results confirm the cur- 
rent standard picture that the universe is dominated by 
dark energy with a cosmological constant-like equation 
of state, independently of type la supcrnovae. 

Although we have quantified our systematics on our 
cosmological results by changing several important as- 
sumptions, there are additional systematic effects that 
could change our quantitative results. One is dust ex- 
tinction by lensing galaxies. Previous studies have shown 
that dust is indeed present in lensing galaxies even if they 
are early-type galaxies, although the measur ed total ex- 
tinct ion is modest, E{B — V) ~ 0.1 mag fe.g.. lFalco et alj 
Il999( l. Since we set the magnitude limit at the i-band 
for which dust extinction is less important than in the 
bluer bands, we expect that this effect will not heavily 
influence our result. Moreover, the flux from the lensing 
galaxy slightly increases the PSF magnitude of the lens 
system, and this effect could compensate the effect of 
dust extinction to some extent. Another effect we have 
neglected is lensing by multiple galaxies. Theoretically 
the probability for such multiple lensing is just a few 
percent of the lensing probab ility by a single galaxy at 
z < 2 (|Moller fc Blainll200lh and therefore can be ne- 
glected. However, it only considered chance superposi- 
tions along the line of sight and ignored nearby correlated 
galaxies that could dominate the contributi on to multi- 
ple lens events (e.g.. lCohn fc Kochaneldl2004f l. Moreover, 
the fact that one o f our lens sample, SDSS JlOOl+5027 
(jOguri et al.ll2005al ). has a secondary galaxy near the pri- 
mary lensing galaxy suggests that the effect needs to be 
addressed carefully. 

In this paper we have assumed a flat universe. Al- 
though this is a reasonable assumption given that 
virtually all current cosmological const raints are con- 
sistent with a flat universe ( e.g., [T egmark et alj 
2006; Ichikawa fc Takahashi 2 0071: ISperge l ct al. 200'J 
Wang fc Mukheriedl2007l : iAllen et al.ll2007r ). it is stiU of 
interest to consider non-flat universes. In particular the 
introduction of both non- flatness and w ^ ~1 results in 
additional stro ng d egeneracy between cosmological pa- 
rameters (e.g., iLind er 2005), and thus may require an 
additional independent cosmological probe in order to 
obtain tight constraints on individual parameters. 

Cosmological constraints presented in this paper are 
obtained from the SQLS DR3 sample, which represents 
- 40% of the fuU SDSS data. The extension of the SQLS 
lens sample to the SDSS DR5 and SDSS-II is in progress, 
therefore we expect that the statistical errors will im- 
prove significantly in the near future. Better constraints 
may be obtained by considering the SQLS and CLASS 
samples jointly. In addition to the extension of the lens 
sample, it is of great importance to reduce the system- 
atic errors. An important advantage of our optical survey 
over radio surveys is that there is no systematics from the 
source redshift distribution, which was the biggest source 
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of systematic error in the CLASS analysis (|Chael 120031) . 
The velocity function of galaxies and the faint end quasar 
luminosity function are expected to converge in the near 
future as current large-scale surveys are completed, thus 
we expect the systematic errors can be reduced in future 
analyses of the SQLS lenses. If we have a large enough 
sample of lenses, we may be able to reduce the systematic 
effect further by calibrating the velocity function from 
the observed image separation distribution itself. 

Although in this paper we have restricted ourselves to 
galaxy-scale lenses to study dark energy, our statistical 
lens sample contains group- or cluster-scale lenses as well. 
The number of cluster-scale lenses in our sample is quite 
sensitive to the abundance of clusters at intermediate 
redshifts (z ^ 0.5), and therefore was used to study as 
(|Oguri fc Keetonll200i : ILi et al.l[200l . The full image 
separation distribution from galaxy- to cluster-scales will 
be valuable in understanding how galaxies are populated 
in dark matter halos. 
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